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Coulomb Correlations and Pseudo-gap Effects in a Pre-formed Pair Model for the 
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We extend previous work on pre-formed pair models of superconductivity to incorporate Coulomb 
correlation effects. For neutral systems, these models have provided a useful scheme which interpo- 
lates between BCS and Bose Einstein condensation with increasing coupling and thereby describes 
some aspects of pseudo-gap phenomena. However, charge fluctuations (via the plasmon, uo p ) signifi- 
cantly modify the collective modes and therefore the interpolation behavior. We discuss the resulting 
behavior of the pseudo-gap and thermodynamic quantities such as T c , \ an d C v as a function of uj p . 



The role of the pseudo-gapH in the high T c cuprates 
is emerging as an important indicator of the nature of 
the superconductivity as well as the normal state. There 
are two widely discussed but competing explanations for 
pseudo-gap effects but no clear and decisive experiments 
to support one scenario over the other. Early observeu 
tions associated the pseudo-gap with magnetic pairingcl 
above T c (often called the "spin gap"). It is now clear, 
however, that some form of normal state pairing is seen 
in photoemission as well as charge transport data. More- 
over, at least in the photoemission data, the pseudo-gap 
appears to have the d-wave symmetrya of the ordered 
state and this leads naturally to the association of this 
"gap" with precursor superconductivity.H~El This second 
scenario is further supported by the observation of low 
dimensionality and short coherence lengths in high T c su- 
perconductors, which suggests important deviations from 
ideal mean field or BCS transitions. Indeed, the approach 
of the present paper assumes the precursor superconduc- 
tivity scenario, in large part because it is important to 
establish, at least as a base-line, the extent to which such 
superconducting "fluctuation" effects may be responsible 
for pseudo-gap behavior. 

Among those models which subscribe to a precur- 
sor superconductivity scenario there are additionally two 
rather distinct viewpoints. Emery and Kivelsonfl have 
argued that the pseudo-gap state of the cuprates is sim- 
ilar to that observed in granular films where phase co- 
herence is not fully established, although large regions 
of the material have a well established superconducting 
amplitude. Because it is small, in some sense, in the 
cuprates their approach focuses on n/m* or alternatively 



"phase stiffness" 
have focused on 



on the plasma frequency ui p as tl 
parameter. Alternatively, other: 
the observed small size of the superconducting correlation 
length £ to argue for important corrections to BCS tin 
ory associated with pre-formed or nearly-formed pairsl 
which exist well above T c and therefore give rise to sig- 
nificant pseudo-gap effects. The present paper is based 
on the viewpoint that in the cuprates the characteris- 
tic parameter of the charge degrees of freedom, n/m* or 
equivalently oj p , should be treated on a relatively equal 
footing with the correlation length, £. 



To study the role of Coulomb interactions on pseudo- 
gap phenomena, we adopt a natural microscopic frame- 
work which incoporates charge fluctuations into theories 
which treat the cross-over from BCS pairing to Bose- 
Einstein condensation (BEC) of pre-formed pairs. In 
neutral systems, this gross-over has been studied by a va-. 
riety of investigators.Uuu Numerical simulation studies,Q 
which have been performed in the context of attractive 
Hubbard models, include, in principle, all diagrammatic 
contributions. On the other hand, analytical work has 
mostly been confined to the T-matrix approximation. 
The issue of non-conserving and conserving T-matrix 
schemes has been widely discussed in the literatureO in 
the context of the BCS-BEC cross-over problem. In 
the original work of Nozieres and Schmitt-Ritik a non- 
conserving approach was used. Recent workH on neu- 
tral systems has extended this scheme using a T-matrix 
approximation which satisfies global conservation laws 
and in the process introduces renormalized Green's func- 
tions into the generalized susceptibilities. In the charged 
system, as a consequence of gauge invariance, the ana- 
logue renormalized susceptibilities must then appear in 
the particle-hole channel. As has been known for some 
timejj however, the collective mode spectrum is then 
treated incorrectly at this level of approximation and a 
more sophisticated scheme is needed. In order to avoid 
this complexity and to develop an intuitive understand- 
ing of the effects of charge, however, we restrict the anal- 
ysis, in this paper, to the more familiar scheme intro- 
duced by Nozieres and Schmitt-Rink and defer consider- 
ation of a fully conserving formalism. We note, however, 
that our formulation will be locally conserving and in the 
case of charged systems this approximation does not yield 
qualitatively different physics from that expected using 
a globally conserving approach. Furthermore, it is our 
contention that mode-mode coupling effects, which are 
ignored in all T-matrix based schemes, will ultimately 
lead to important insights which we will discuss in a fu- 
ture paper. 

The Hamiltonian under consideration contains an at- 
tractive interaction Vk,k'> parameterized by a coupling 
constant g, as well as long range Coulomb terms. For 
definiteness wc take the same separable pairing poten- 
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tial, 14;, k' = gvkVk' where Uk = (1 + k 2 /k^)^ 1 / 2 ; as was 
used initially by Nozieres and Schmitt-Rink. Within this 
model, the pairing energy scale (or "Debye frequency") 
is the Fermi energy. We assume a three dimensional free 
electron model for the electrons and defer discussion of 
anisotropy effects until later in the text. It is assumed 
that in the cuprates there is sufficient inter-layer hop- 
ping so that a strictly two dimensional model and its 
associated Kosterlitz-Thouless transition is not the ap- 
propriate starting point. We generalize the path integral 
formulation of Ref . [|, replacing the usual fermionic fields 
by Nambu spinors, = (cL-, c_kj.) ; and decoupling in 
four real fields given by r\i = V riV* { T ^^ i = . . . 3, are the 
identity and three Pauli matrices). The interaction in the 
off-diagonal channels (i = 1,2) is the pairing interaction 
while the Coulomb term appears in the diagonal chan- 
nels (i = 0,3). Finally, the field r/o may be eliminated 
by a suitable gauge transformation and the thermody- 
namic potential f2 is computed for the remaining degrees 
of freedom at the Gaussian approximation level, 

tt = f7o + | £ m det (F q i J - ^ ]>> q ( 3 \ 

1,Vrn q 
^ = Sij ~TV® tr fr Sk+q/2,*, T i &-qA<-,-*.» 

Here fio is the usual contribution from non-interacting 
fermions and n t denotes the corresponding number den- 
sity of free fermions. In the interaction = g, v^' = 
w k for % = 1,2, and V^ 3) = Ane 2 /^ 2 , = 1. Finally, 
Gq.oj is the single particle Green's function in Nambu 
space and T is the temperature, while v m and cj; are the 
even and odd Matsubara frequencies, respectively. Above 
T c the fluctuation propagator T is diagonal; thus all three 
channels (the particle-particle, particle-hole and uncorre- 
cted fermions) contribute additively to as well as to 
the various derived thermodynamic quantities which we 
calculate below. n 

In the present formalism,!^ the transition temperature 
is obtained from the BCS gap equation with a self consis- 
tently determined chemical potential, /x, obtained from 
the condition n to t = —dCl/dfi. The resulting coupled 
equations are solved numerically for \i and T c . In the 
limit of arbitrarily large g the effects of the Coulomb 
interaction drop out, since mode-mode coupling is ne- 
glected at the Gaussian level, so that T c is given by 
the ideal BEC temperature. At small g, the collective 
mode contribution to /x becomes arbitrarily small and the 
BCS limit is approached, albeit with a Coulomb renor- 
malized chemical potential. t£l When the superconducting 
coupling constant g vanishes, the above form for f2 re- 
flects the plasmon contribution and reduces to that of 
the usual Coulomb gas.E2l The Gaussian approximation 
to gives an RPA-like treatment of the collective modes 
and so provides a reasonable scheme for interpolating be- 
tween these two limits. It should be noted, however, 
that Coulomb pseudo-potential effects (which would act 



to renormalize g) are, for simplicity, not included in our 
calculations. Here we focus principally on the effects in- 
troduced by the long range Coulomb interaction (which 
enters via the parameter u p 2 = Anne 2 / (m* k)) . 
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FIG. 1. (a) The variation of T c as a function of coupling 
constant in the BCS and Gaussian approximation theories for 
neutral and charged systems with increasing plasma frequen- 
cies, (b) Tc normalized by the Coulomb mean field result (see 
text) To for the same pfasma frequencies as in (a). The inset 
in (b) plots the number of paired eiectrons vs g/g c for the 
same parameter set. 

To illustrate the effects of charge fluctuations we plot, 
in Fig. |l|a, T c as a function oig/g c for various values of the 
plasma frequency. The dotted line represents the BCS 
result (for the neutral system) and the solid line is the 
corresponding (neutral) transition temperature which in- 
cludes Gaussian fluctuations. The remaining curves from 
left to right demonstrate the effects of increasingly large 
Lo p . For the purposes of focussing on charging effects only 
we fix Ep and kp/ka; in this way all the Gaussian derived 
curves have the same high g asymptote. Moreover, with 
these assumptions the neutral system reference curve is 
unchanged as lu p is varied. Two effects of Coulomb in- 
teractions can be observed in Fig. |l|a: (i) T c decreases 
with increasing u> p , and (ii) the non-monotonic behavior 
as a function of g found for the neutral case (which is be- 
lieved to be unphysical)El progressively disappears with 
increasing u> p . 

The first observation, which is perhaps the more im- 
portant, is a consequence of the fact that an attraction 
due to Coulomb interactions in the particle-hole channel 
reduces the effectiveness of the attraction in the pairing 
channel. We illustrate this first point more directly in the 
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inset of Fig. |lp, which plots the number of superconduct- 
ing pairs as a function of increasing uj p . This effect may 
seem counter to the expectation that systems with larger 
uj p will have reduced T c suppression (i.e. larger T c ) due to 
phase fluctuations. However, when the appropriate refer- 
ence temperature is used, Coulomb interactions are found 
to lead to better agreement with mean field theory; in 
this sense superconducting fluctuations are, indeed, sup- 
pressed by Coulomb interactions. We plot in the main 
portion of Fig. [l|b the ratio of T c to the critical temper- 
ature Tq obtained by neglecting pair fluctuations (but 
including Coulomb effects) for the same range of plasma 
frequencies as in Fig. |l|a. As can be seen, the larger the 
plasma frequency, the higher the ratio T c /Tq and thus 
the better the agreement with a mean field treatment of 
the pairing channel. 



a) CjT vs T 




T 



FIG. 2. Temperature dependent specific heat (a) and spin 
susceptibility (b) for the same parameters as in Figures 1. 
The dotted lines are the BCS results (with slightly tempera- 
ture dependent chemical potential) . The inset in (a) plots the 
total number of pairs vs T, along with the number at q = 0. 

A more convenient way of illustrating fluctuation ef- 
fects and the associated role of the plasma frequency, 
however, is to study thermodynamic properties directly. 
Here pseudo-gap effects enter as precursor superconduct- 
ing contributions in, for example, the specific heat and 
spin susceptibility. As in Fig. [l|b, these fluctuation effects 
are expected to weaken as the plasma frequency increases 
and mean field behavior is restored. We plot C v and \ 
in Fig. |2| for the same parameter set as in Fig. [I], with 
g/g c set equal to unity. This choice of coupling strength 
is consistent with the observation of relatively short co- 
herence lengths in the cuprates and with the claim that 



high T c compounds lie close to the bound state limit. □ 
Here the dotted lines again represent the BCS value 
(with a properly T dependent fi). As is appropriated 
for Gaussian fluctuation theories, the specific heat varies 
as (T — Tc)" 1 / 2 . Comparison with the neutral (solid) 
reference curve shows that the effects of variable uj p are 
not as evident in the specific heat as in the spin suscep- 
tibility. Similarly, deviations from mean-field behaviour 
appear at higher temperatures in \ than in C v . 

It is worth noting that in the present formulation there 
is no well defined "onset temperature" for the appearence 
of a pseudo-gap, as might exist if a sharp phase transition 
occured at some temperature T* above T c . Moreover in 
discussing the onset of fluctuations, it should be stressed 
that their appearance is unrelated to the temperature 
dependence of the total number of fluctuations or pre- 
formed pairs. As shown in the inset of Fig. ||a, for fixed 
w p , the total number of pairs is relatively constant up to 
temperatures many times higher than T c ; however, their 
distribution shifts to lower momenta as T c is approached 
and long range coherence is established. Detailed calcu- 
lations indicate that varying lo p does not alter the critical 
behaviour in the pairing channel as T c is approached; the 
narrowing of the pseudo-gap region with increasing uj p is 
a result of the smaller relative contribution that pairing 
fluctuations make to the thermodynamics as Coulomb 
correlations become more dominant. Finally, it should 
also be noted that the neutral system yields unphysical 
thermodynamic behavior at smaller coupling constants 
than when Coulomb effects are included. This is illus- 
trated by the negative values of \ indicated in the figure. 
Such unphysical behaviour has been shown to result from 
the non-conserving nature of the approximations used.t2l 

While the above figures were designed to illustrate the 
effect of Coulomb correlations, they do not fully rep- 
resent the physical system in which variations in the 
plasma frequency are necessarily associated with changes 
in the electronic energy scale. In reality both uj v and 
Ep = kp 2 /2m* depend on similar combinations of the 
carrier density n and effective mass m*. As the insu- 
lator is approached uj p decreases as lo p 2 ~ x (where x 
denotes the number of doped holes); however, whether 
one assumes a Fermi liquid (n » 1 + i) or non-Fermi 
liquid (n ss x) approach to the insulator, it follows nec- 
essarily that the electronic energy scale Ep must also 
vanish as the hole concentration x approaches 0. In sce- 
narios based on electronic pairing mechanisms, therefore, 
it is difficult to escape the conclusion that the onset tem- 
perature for coherent superconducting fluctuations, T*, 
should also become small as the insulator is approached. 
Our numerical calculations of C v and x exhibit this effect, 
principally because ours is a single energy scale theory: 
both T* and T c derive from the same pairing mechanism. 
This behavior is in contrast to experimental where even 
in highly underdoped systems T* is of the order of 100K 
or more. Large T* seems to be most naturally associated 
with a high energy scale in the insulating parent com- 
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pound, such as a magnetic energy.EI However, in such a 
scenario it then becomes problematic to understand how 
the other important energy scaleH to p enters to determine 
T c . We speculate that a cross-over from three to two 
dimensionality may play some role in lowering T c at low 
doping concentrations where enhanced-(quasi-2d) critical 
fluctuation effects are most apparent O 
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FIG. 3. Variation of T c with lo v (in units of eV, where 
optimal doping corresponds to uj p ~ 1.2eV) with Ep held 
fixed (curves (a) and (b)) as well as using self-consistently 
determined Ef (curve(c)). Curve (b) corresponds to 2 + e 
dimensions with e = when ui p = and e = 1 near optimal. 
See text for details. 

In Fig. [| we explore these issues within the context 
of our model. We consider three situations in which the 
Gaussian derived T c is plotted as a function of to p (which 
may be directly related to hole concentration x) to arrive 
at a form of "phase diagram" . The solid line corresponds 
to the case in which the characteristic energy for pairing 
(called il Ep") is held fixed, and as observed in Fig. 0a, 
T c monotonically decreases. This should be contrasted 
with the situation in which Ep is allowed to vary self- 
consistently (dotted line) in accord with the measured 
ujp. Here for dehnitencss, we assume^ Luttinger volume 
Fermi surface. The latter case yieldsa T c « oj p \ however, 
T c is suppressed at low doping primarily as a result of the 
lowering of the electronic energy scale, Ep, rather than 
from increased phase fluctuations. Neither of the above 
cases is entirely satisfactory for understanding the larger 
pseudo-gap regime at low doping. We address this issue 
phcnomenologically by introducing the effects of a dimen- 
sionality crossover for fixed pairing energy scale Ep and 
calculating T c for a system in 2 + e dimensions (dashed 
curve) where e varies smoothly from zero when lo p ss 
(at half filling) to one when io v w 1.2eV (corresponding 
to optimal doping). Here T c is suppressed at low dop- 
ing, as a consequence of 2d fluctuation effects. Moreover, 
one can associate T* with the solid curve, which exhibits 
the observed experimental trends. Thus the low doping 
regime is characterized by a large pseudo-gap. On the 
other hand, at higher doping or lo p , both T c and T* con- 
verge and the pseudo-gap region vanishes. Whether or 
not this phenomenology is appropriate, the above discus- 



sion underlines the importance of multiple energy scales 
(w p , as distinct from the pairing energy scale) and the 
possible role of a dimensionality crossovert2l in under- 
standing pseudo-gap behavior. 

In summary, we have presented a pre-formed pair 
model in which the effects of Coulomb correlations are 
clearly seen to suppress superconducting fluctuations 
(above T L ) and thereby tune pseudo-gap behavior in the 
calculated specific heat and spin susceptibility. More- 
over, these two thermodynamic variables, as a function 
of T, are found to be reasonably consistent with exper- 
iment. However, it should be stressed that within our 
microscopic model, the effects of Coulomb correlations 
enter in a rather different way than has., been assumed 
in previous phenomenological schemes. ltEj Our approach 
focusses more directly on correlated pairs rather than 
superconducting grains; therefore, phase and amplitude 
fluctuations appear on a relatively equivalent basis. As a 
consequence, introducing Coulomb interactions into the 
pre-formed pair formalism leads to a narrowing of the 
pseudo-gap region by providing a competing attraction 
in the particle-hole channel and thus reducing the effec- 
tiveness of superconducting pairing. 
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